Direct and inverse pumping in flows with homogeneous and non-homogeneous swirl 
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The conditions in which meridional recirculations appear in swirling flows above a fixed wall are 
analysed. In the classical Bodewadt problem, where the swirl tends towards a fixed value away 
from the wall, the well-known "tea-cup effect" drives a flow away from the plate at the centre of the 
vortex. Simple dimensional arguments applied to a single vortex show that if the intensity of the swirl 
decreases away from the wall, the sense of the recirculation can be inverted, and that the associated 
flow rate scales with the swirl gradient. Only if the flow is quasi-2D, does the classical tea-cup effect 
take place. This basic theory is confirmed by numerical simulations of a square array of steady, 
electrically driven vortices. Experiments in the turbulent regimes of the same configuration reveal 
that these mechanisms are active in the average flow and in its fluctuating part. The mechanisms 
singled out in this letter provide an explanation for previously observed phenomena in electrolyte 
flows. They also put forward a possible mechanism for the generation of helicity in flows close to 
two-dimensionality, which plays a key role in the transition between 2D and 3D turbulence. 



INTRODUCTION 



The teacup effect is one of the mechanisms actively 
mixing sugar in a cup of tea when stirring it with a 
spoon. Under the effect of rotation, a centripetal pressure 
gradient builds up in the fluid to oppose the centrifugal 
force. In the boundary layer near the bottom of the cup, 
the flow is slow so the centrifugal force collapses and the 
same pressure gradient drives a convergent flow towards 
the centre of the cup. This feeds the meridional recir- 
culation, which is actually responsible for mixing (this 
configuration is that of the Bodewadt problem). The 
configuration where a solid wall rotates under a still fluid 
(Ekman problem) leads to a reversed meridional flow for 
the same reasons [l|, [2[ . 

This type of mechanism is at play in a number of less 
anecdotal processes, both industrial and natural, such as 
the stirring of liquid metals (in the Bridgeman process 
to grow silicone crystals for example [3j|) or the genera- 
tion of cyclones, where it controls the redistribution of 
momentum, heat or chemicals in the flow. It is particu- 
larly important in experiments aiming to reproduce two- 
dimensional turbulence or to understand the transition 
between two and three-dimensional turbulence [J, |5| : the 
authors of Refs. [g, [7J recently discovered that the ap- 
pearance of a third velocity component played a role in 
the break-down of the inverse energy cascade of 2D tur- 
bulence, and was therefore central to understanding this 
transition. Recent numerical simulations |8j were able to 
link the existence of the inverse cascade to the presence 
of helicity even in 3D turbulence. Since the teacup effect 
is precisely a source of helicity linked to the presence of a 
bottom wall, it most likely plays an important role in the 
transition between 2D and 3D turbulence dynamics. The 
authors of Ref. [9(| attempted to suppress it by inserting 
a "buffer" fluid layer between the container wall and the 
layer of fluid where turbulence was forced, only to dis- 



cover that secondary flows still subsisted. Their presence 
was attributed to confinement itself but also to the non- 
homogeneity of the forcing, which consisted of passing 
an electric current through the fluid layer (a conductive 
electrolyte) , placed over an array of magnets of alternate 
polarity. Remarkably, pumping in this configuration was 
inverse, with the fluid diving to the centre of vortices in 
the core of the flow, whereas in homogeneously forced 
flows between two planes, direct pumping is expected as 
in the teacup effect [lOj. Thus, although the teacup ef- 
fect is well understood, the conditions in which secondary 
flows appear and simply which way they flow is not yet 
clear. 

In this letter we put forward a mechanism to explain 
how either direct or inverse pumping arises, depending 
on the homogeneity of the forcing. This scenario is tested 
against numerical simulations and experiments on a flow 
of liquid metal between two parallel planes, subject to 
a transverse magnetic field Be, and where the flow is 
driven by injecting electric current at one of the walls. 
This setup offers a convenient way to control the homo- 
geneity of the forcing, since for strong magnetic fields, the 
Lorentz force diffuses momentum across the fluid layer 
over a time scale of (H/l±) 2 p/(crB 2 ), thus making the 
flow quasi-two dimensional, and homogeneously forced 
across the fluid layer (H, l±, p and a are the thickness 
of the fluid layer, the size of the swirling structure, the 
density and electric conductivity of the fluid). Direct 
pumping was shown to occur in this case 1 1 lM 1 3f ! At 
more moderate magnetic fields, momentum diffusion can 
be made slower than eddy turnover time l±/u(l±_) and 
the forcing intensity decreases away from the wall where 
electric current is injected, thus becoming inhomogeneous 
[14( (u(l±) is a typical velocity within the structure). 



GOVERNING EQUATIONS 

The velocity and pressure fields u and p of an incom- 
pressible flow of a fluid of density p, viscosity v are gov- 
erned by the Navier-Stokes and continuity equations: 



(d t +u- V)u + -Vp = iA7 2 u - 
P 
V-u = 0. 



(1) 
(2) 



To illustrate the phenomena of direct and inverse pump- 
ing, we shall rely on the example of electrically forced 
flows. The principle is to apply an externally generated, 
homogeneous magnetic field Be z to an electrically con- 
ducting fluid and to inject electric current in one or sev- 
eral points of an otherwise electrically insulating wall or- 
thogonal to e z (here at z = 0). At each such electrode, 
the Lorentz force creates a vortex spinning around e z and 
extends all the further in the core as B is high |l4 , llol | . 
Assuming the magnetic Reynolds number Rm = paUL 
remains small (U, L are typical velocities and length, p 
is the magnetic permeability of vacuum and a the elec- 
tric conductivity of the fluid.), then the magnetic field 
associated to the current J within the fluid is negligible 
compared to B, and the Lorentz force to which the flow 
is subjected expresses as f = Bp~ l 3 x e z . J is coupled 
to u through Ohm's law and charge conservation: 



J = a(— V(f> + Bu x e 2 
V- J = 0, 



(3) 

(4) 



where <j> is the electric potential. The generic geome- 
try is that of a channel with no-slip, impermeable and 
electrically insulating walls located at z = and z = H. 
The flow is governed by two non-dimensional parameters: 
the Hartmann number Ha = BH\J a /(pv) and and the 
Reynolds number Re = Y/v. Y = I / '(2ir ^j 'a pv) is the 
circulation that would be induced around an electrode, 
in the plane just outside the Hartmann boundary layer at 
z = 0, by injecting a DC current of intensity / through 
it, if no viscous dissipation was present |l4 , Il5| . Ha 2 
represents the ratio of Lorentz to viscous forces but Ha 
and Re can be thought of as non-dimensional measures 
of B and / respectively. As mentioned in introduction, 
Ha controls the time scale of momentum diffusion along 
B and so decreasing it incurs steeper velocity gradients 
along e 2 in the core of the flow, and increases the inho- 
mogeneity of the swirl. 



NUMERICAL SYSTEM 

In the numerical simulations, the fluid is confined in 
a rectangular box of size L 2 x H with L = 0.06 m 
and H = 0.1 m. The working fluid is GalnSn, an eu- 
tectic alloy of Gallium Indium and Tin, of conductivity 
a = 3.6xl0 6 fi -1 .m -1 , viscosity i/ = 4xl0 -7 m 2 .s _1 and 



density p = 6.4 x 10 3 kg.m -3 . The container is placed 
in a uniform magnetic field Be z so that besides the two 
Hartmann walls orthogonal to it, there are also four walls 
parallel to it. All walls are impermeable, no slip and elec- 
trically insulating. The flow is driven by injecting a DC 
current or alternate polarity through a square array of 6 
x 6 electrodes of diameter 1 mm, mounted flush at the 
bottom wall, spaced by distance Li = 0.01 cm. 
The code is based on the finite volumes method imple- 
mented in the OpenFOAM framework and solves the 
three-dimensional, time-dependent equations in a segre- 
gated way. The numerical scheme is the consistent and 
conservative algorithm put forward in Ref. [16(. The 
code is described and fully tested in the more complex 
configuration of the flow around a 3D obstacle in Ref. 
[17| . To summarise it, the spatial discretisation is of sec- 
ond order and the time-scheme is a second order implicit 
pressure-velocity formulation and the pressure-velocity 
coupling is solved with the PISO algorithm implemented 
as in Ref. [18j. The time step is chosen so that the 
maximum Courant number C = U At/ Ax is well below 
1 and the maximum of D = vAt/ Ax 2 remains below 
10 (Courant-Friedrich-Lewy conditions). Collocated and 
structured meshes made of 1 917 971 and 4 227 768 rect- 
angular elements are respectively used for simulations at 
Ha = 800 and Ha = 1822. In order to keep the mesh 
orthogonal, the electrodes, which are circular in the ex- 
periment, are modelled by squares of size 1 mm. The 
mesh is regular outside of the boundary layers and re- 
fined near the walls so as to always keep 4 points across 
the Hartmann layers and 3 points across the Shercliff lay- 
ers. 

For each value of Ha, the flow is initially at rest and the 
lowest current is injected. In all computations presented 
here, the flow stabilised into a steady state, which was 
used as initial condition for the simulation at the same 
Ha, for the next value of Re up. 



MECHANISMS OF INVERSE AND DIRECT 
PUMPING 

We shall first illustrate the basic mechanism of inverse 
pumping on the example of a steady flow in the con- 
figuration described in the previous section. Figure Q] 
shows the pressure contours and streamlines obtained at 
Ha = 1.822 x 10 3 and Re = 171. These are typical 
of the flow patterns obtained for other parameters. The 
flow consists of a square lattice of 6 x 6 vortices rotating 
along e 2 with alternate spin. A strong downward flow 
exists at the centre of each of them that loops back up 
in the outer part of the vortices. To illustrate the under- 
lying mechanism, we shall reason on a single vortex of 
radius l±, with associated polar coordinates centred on 
it (r 2 = (x + Li/2) 2 + (y - Li/2) 2 , and the radial and 
azimuthal velocity in this vortex respectively correspond 
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FIG. 1. Streamlines ol (u x ,u z ) in the plane y = Li/2, in- 
tercepting the axis of rotation of two vortices and contours 
of pressure (colours). Pressure is normalised by pu 2 M where 
um = m&x{u y (x,z)} is the maximum azimuthal velocity 
within any given vortex in the vessel. Downward vertical jets 
occur right in the region of lowest pressure at the centre of 
vortices at x = ±Li/2. 



to u x and u y in the plane y = Li/2, on figure [l}. The 
pressure contours in figure[T]and vertical profiles of p and 
u z (figure [21 left) show that the downward jet coincides 
with a strong pressure gradient. Plots of the different 
terms in ([T])-e z (figure [2], right) show that the latter is 
only opposed by viscous friction: 



d z p 



Uz 



(- r >) 



The pressure drop at the centre of the vortex results from 
the centripetal pressure gradient that opposes the cen- 
trifugal force. 

p(r = 0, z) - p(r = l±,z) ~ -pu 2 e {r = l±,z). (6) 

Contours of p on figure [1] suggest that the pressure gra- 
dient outside of the vortex d z p(l±,z) is negligible. The 
maximum u^f of u z (Q,z) is located at z = Zq, near the 
bottom wall, but still outside the boundary layer (see 
figure [2 left) and can be estimated from ([5]) and ([6]): 



u z (0,z o ) 



i?e N 



(7) 



ue(l±,zo) 

where i?e v = 2d z ue(l±, zq)1\/v is a Reynolds number 
built on the vertical gradient of the swirl. Eq. ([7]) ex- 



-p(-l/2,z) 
■u z (-l/2,z) 
■•Mj, (-0.372, 2) 
'"0* (-0.372, 5) 
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FIG. 2. Left: profiles of pressure and velocity in the 
plane y/Li = 1/2, along the columnar vortex centred at 
(x/Li,y/ Li) = ( — 1/2,0). Right: variations of the different 
terms in Jl])-e z along the axis of rotation of the same vortex. 
The vertical pressure gradient balances lateral viscous friction 
to a good approximation in the core. Velocities, pressure, x 
and z coordinates marked q were respectively normalised by 
um = max.{u y (x, z)}, pu Ml Li and H. 



presses that vertical motion is driven by the inhomogene- 
ity of the swirl. Mass conservation requires that the flow 
impacting the wall must turn radially: 



u r (z ) 



— — U z (0, Zn) 

2h r v ' 



-u e (lj_,z )—Re v , 



(8) 



where h r is the height over which the return flow takes 
place, typically of the same order as l± in all simulations. 
Unlike inverse pumping, direct pumping is driven in the 
wall boundary layer, by the radial pressure gradient that 
builds up in the core to oppose the centrifugal force ([6]). 
If S is the boundary layer thickness, ([I])-e z implies that 
the pressure cannot vary across it and so the radial pres- 
sure gradient there is the same as in the core. Since ug 
becomes small in the boundary layer, so does the centrifu- 
gal force pug/r. Hence, (JTJ-e,- there expresses that in the 
boundary layer, the centripetal pressure gradient is bal- 
anced by viscous friction. If inverse pumping is present, 
however, the associated radial flow in the core, expressed 
by ([5]), considerably enhances radial friction vd^ z u r in the 
boundary layer. Denoting quantities within the bound- 
ary layer with a superscript b, it comes that 



d rP b 



pv_ 



(u b r -U r (zo)), 



(9) 



and from © and ©, 



Ufl 



2h r 



Re^ 



Re, 



(10) 



where, unlike Re°, the Reynolds number Re = 
ue{zo)l_\_/i> is based on ug(zo), a measured quantity. For 
direct pumping to exist, the flow must be directed radi- 
ally inwards within the boundary layer. When d z ug < 0, 
as in the case studied here, this imposes the condition: 



if) r 



(id 



In all numerical simulations, h r was typically of the order 
of l±. Boundary layers are thin, so l±/8 >> 1 (in electri- 
cally driven flows the wall boundary layer is a Hartmann 
layer, with thickness 8/H = Ha -1 , so l±/8 ~ Ha). Since 
the flow is quasi-2D as soon as d z u = 0(8/l±) outside 
the boundary layer, (fTTj) can only be satisfied when the 
swirl is quasi-2D. This result explains why direct pump- 
ing was never observed in any of our simulations: in all 
numerically accessible regimes, the inhomogeneity of the 
swirl was always sufficient to drive inverse pumping. Sim- 
ilarly, in the wall bounded electrolytes layers studied by 
[19j . vortices were driven by imposing an electric cur- 
rent in the field created by permanent magnets placed 
underneath the layer. The swirl inhomogeneity resulted 
directly from that of the magnetic field and so d z ug and 
ug/H were of the same order, making it impossible for 
(TTTJ) to be satisfied. It is therefore not surprising that 
inverse pumping was observed in this case too. 
In quasi-2D flows, by contrast, Re v ~ and from (p~0|) 
and (|2|). direct pumping drives a flow from the boundary 
layer into the core: 



Ufl 



Re. 



(12) 



In electrically driven flows, 8 = H/ Ha and sc aling ([T2 



coincides with the expression derived by Ref. |ll| from 
matched asymptotics. A generalisation of this result 
was recently derived that includes both magnetohydro- 
dynamic flows and flows with a background rotation [13| . 
It should also be noticed that should the swirl be inho- 
mogeneous in such a way that <9 2 |«e| > 0, then direct 
pumping would occur, but u z would still scale as in ([7]). 



INVERSE PUMPING IN THE SQUARE VORTEX 
ARRAY 



Numerical simulations of the vortex array present spe- 
cific features that differ from those of an isolated vortex, 
on which the theory is based. Despite this difference, 
figure [3] shows that the vertical velocity u z (zq) within 
a single vortex in the array linearly increases with i?e v . 



Scaling (JT]) , is thus still satisfied, albeit with a reasonable 
amount of data scattering. This scattering can be partly 
attributed to the plotted values being local ones, and 
partly to interaction between vortices. Most importantly, 
data obtained for different values of Ha collapse on the 
same line, which confirms i?e v as the single relevant pa- 
rameter for inverse pumping. In other words, it is the 
gradient of swirl that determines inverse pumping, not 
how this gradient is generated, electromagnetically or in 
any other way. Furthermore, all simulated cases involve 
strongly inhomogeneous swirl, with <9 2 |u£»| > \ug\/H and 
exhibit inverse pumping, in agreement with condition 
(fTTj) too. 

Unlike when isolated, vortices embedded in an array un- 
dergo a strong influence from neighbouring vortices. In 
an infinitely extended square array, this would mainly 
translate into a loss of axisymmetry of individual vor- 
tices, with streamlines progressively deforming to a 
square shape, away from the vortex cores. In our numeri- 
cal simulations, by contrast, the 6x6 array is bounded by 
lateral walls. These incur significant friction on periph- 
eral structures, which are consequently weaker, and less 
influent than those near the centre. As a consequence, 
the upper part of the vortex axes is slightly diverted away 
from the centre, instead of being straight. Nevertheless, 
as vortices are driven by injecting current from the bot- 
tom wall, they remain attached to the electrodes in this 
region, and the effect of the imbalance between centre 
and peripheral vortices becomes more visible away from 
the bottom wall (This can be noticed in figure [1]). A 
second property of vortex arrays is that radial flows in- 
duced in neighbouring vortices by inverse pumping col- 
lide to form a strong return flow that spirals up in the 
region between vortices (This pattern is reflected in the 
vertical jet at x = in the streamlines of (u x , u z ) in fig- 
ureQ]). When Re is increased, these phenomena become 
more pronounced, and at Re = 171 (for Ha = 800) and 
Re = 512 (for Ha = 1822), recirculations merge in the 
upper part of the vessel, to the point where they cannot 
be distinguished from each other anymore. 



DIRECT AND INVERSE PUMPING IN 
TURBULENT FLOWS 

We shall now illustrate the occurrence of inverse pump- 
ing in more complex flows, such as turbulent flows. Since 
the corresponding regimes lie beyond the reach of nu- 
merical simulations, these shall be analysed experimen- 
tally. The experimental setup closely matches the config- 
uration of the numerical simulations with one difference: 
the dimensions of the vessel across the magnetic field are 
L x L = O.lmxO.lm instead of L x L = 0.06mx0.06m. 
The distance between electrodes where the current is in- 
jected can be set to Li = 0.01m (as in the numerical sim- 
ulations) or Li = 0.03m. A full description of the rig can 
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FIG. 3. Scaled maximum vertical velocity vs^ _Re v , from 
numerical simulations of the 6x6 vortex array in steady state 
(Re v gives a measure of the gradient of swirl along e z ). Re 7 
was calculated with l± = Li/4. 



be found in [1J, |20j. The influence of the walls becomes 
more pronounced when Li/ L increases. In this regard, 
the numerical simulations present an intermediate case 
between the two injections scales available in the experi- 
ment. For the purpose of this work, the rig was equipped 
with ultrasound sensor-transducers fitted flush in the top 
wall and in one of the lateral walls, connected either 
to a DOP1000 or a DOP3010 ultrasound velocimeter, 
manufactured by SIGNAL PROCESSING. Ultrasound 
velocimetry is the method of choice to obtain instanta- 
neous velocity profiles in opaque fluids and is now well 
developed for liquid metal flows |2l| 
provide instantaneous profiles of u z 
(—1.5cm, 1.5cm) (probe VI, aligned with a current injec 



22] . These probes 
along e 2 at (x, y) = 



tion electrode) , as well as 5 horizontal profiles of u x along 
e x at z G {0.12H, 0.3 IF, 0.5H, 0.69H, 0.88H} 
(probes HI to H5), also halfway between two electrodes. 
The signals were reliable up to resolutions of 15Hz and 
2mm. In subsequent experiments, a constant current is 
injected in a fluid initially at rest and all results pre- 
sented thereafter are obtained when the flow has reached 
a statistically steady state. The flow is in a turbulent 
state where fluctuations exceed the intensity of the aver- 
age flow. 

Figure |4] shows the time-averaged profiles of 
u z (— 1.5cm, 1.5cm, z) obtained from probe VI, and 
u x (0,0,z) from probes HI to H5, as well as the RMS of 
the fluctuations of these quantities (respectively denoted 
as (•) and (•' 2 ) 1 / 2 ). For L t /H = 0.1, Ha = 1.0932 x 10 4 
and Re = 2.488 x 10 4 , the average and the fluctua- 
tions of the swirl are clearly 3D, as |u K (z)| noticeably 
decreases away from the bottom wall. Inverse pumping 
is present in the lower half of the vessel, in agreement 
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FIG. 4. Profiles of time-averaged horizontal and vertical 
velocities (left) and RMS of velocity fluctuations (right) 
obtained by ultrasound velocimetry. Average velocities 
( resp. fluctuations) are normalised by max{{u x (x, z))}, 
( resp. m&x{(u x (x,z)' 2 ) 1 ' 2 }). Dash: vertical velocity, dia- 
monds: value of u x (average and fluctuations) measured at 
(x/H,y/H,z/H) — (0,0, z/H), dash-dot: order 3 polyno- 
mial fit for Ux, top: Ha = 1.0932 x 10 4 , Re = 2.488 x 10 4 , 
Li/H = 0.1, bottom: Ha = 7.288 x fO 3 , Re = 3.11 x 10 3 , 
L z /H = 0.3. 



with the prediction of ([TTj) . The general shape of profile 
(u z (z)} qualitatively follows 0, in the sense that (u z (z)) 
is more intense where the gradient of swirl is stronger. 
Furthermore, faint direct pumping can be noticed in the 
vicinity of the upper plate (with the flow directed away 
from the plate) . This local effect takes place in a region 
where both the swirl and its vertical gradient are weak. 
Over a short range of values of z, between the regions 
of direct and inverse pumping, u z (z) becomes slightly 
positive, suggesting that a counter-rotating "connecting" 
recirculation is present between the co-rotating recircu- 



lations associated to direct and inverse pumping. The 
profile of fluctuations does not provide the direction of 
the vertical flow, but still shows that these are strong. 
Their intensity is stronger in regions of stronger vertical 
gradients of (m x (0,0, z)' 2 ) 1 ^ 2 , which suggests that they 
are also driven by the inverse pumping mechanism. 
For Li/H = 0.3, Ha = 7.288 x f 3 and Re = 3.11 x 10 3 , 
larger, slower vortices are more sensitive to momentum 
diffusion by the Lorentz force and the flow is corre- 
spondingly closer to quasi two-dimensionality. The av- 
erage flow nevertheless still displays the trace of the 
forcing, with velocities that are highest near the bot- 
tom wall, high near the top wall (because of the strong 
electric current present in the top Hartmann layer) and 
weaker in the middle, where inertial effects take away en- 
ergy of the main flow (The mechanisms governing three- 
dimensionality in wall bounded MHD flows are anal- 
ysed in detail in [14J). Remarkably, this small three- 
dimensionality is sufficient to drive inverse pumping in 
the vicinity of both top an bottom walls and supersede 
direct Bodewadt pumping, in line with our theoretical 
prediction that direct pumping can only occur in quasi- 
2D flows. The profile of (u x (0,0, z)' 2 ) 1 / 2 is practically 
quadratic and symmetric with respect to z/H = 1/2. 
The maximum intensity of the fluctuations corresponds 
to the minimum of the average flow, from which they 
draw energy (near z/H = 1/2). This phenomenon, called 
"barrel effect" was first predicted theoretically [11|, found 
numerically [23j and was recently shown to act as a gen- 
eral mechanism of appearance of three-dimensionality in 
wall bounded flows [13(. These measurements constitute 
the first experimental observation of this effect. Again, 
he vertical velocity is strongly fluctuating, but the low 
value of the gradients of {u x (0, 0, z)' 2 ) 1 / 2 does not allow 
us to conclude as to whether the fluctuations correspond 
to direct or inverse pumping in this case. 



CONCLUSION 

Theory, numerical simulations and experiments concur 
to show that direct recirculations can be inverted (or 
reinforced) when velocity gradients appear in the third 
direction. This phenomenon does not only take place in 
steady vortices but also in fluctuating structures, which 
raises the question of its relevance to fluctuations at the 
different scales of turbulent flows. The appearance or 
the inversion of secondary flows in turbulent structures 
indeed provide mechanisms to create helicity or to 
reverse its sign. Since the sense of the energy cascade is 
tightly linked to signed helicity [8J], they could play an 
important role the transition between direct and inverse 
energy cascades. Nevertheless, our early experimental 
observations of turbulent flows show that inverse pump- 
ing combines with other effects to drive complex flow 



patterns with an helicity that can change sign along 
the transverse direction in a far from straightforward way. 
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